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Abstract 

The  physical  and  electronic  structure  of  the  (1  0  0),  (0  1  0),  (0  0  1)  and  (1  0  I)  faces  of  P-Ga203  are  addressed  using  ah  initio  theory. 
Restricted  Hartree-Fock  calculations,  with  large-core  Ga  and  O  pseudopotentials,  are  done  to  optimize  the  structure  of  first  the  hulk  and 
then  of  slabs  “cut”  in  the  required  orientations.  The  slab  unit  cells  are  fully  relaxed  during  optimization,  and  the  displacements  of  all 
atoms  from  the  ideally-terminated  positions  are  obtained  as  functions  of  depth  into  the  bulk.  For  the  relaxed  slabs,  single-point  density 
functional  theory  calculations  using  the  B3LYP  functional  and  all-electron  basis  sets  are  performed  to  obtain  surface  energies,  ionic 
charges  and  bond  overlap  populations.  All  surfaces  exhibit  a  decrease  in  surface  energy  upon  relaxation,  and  the  local  bonding  at 
the  surface  is  analyzed  by  comparing  nearest-neighbor  bond  lengths  and  overlap  populations  with  those  in  the  bulk.  The  (1  0  I)  surface, 
which  exhibits  a  high  energy  when  ideally  terminated,  undergoes  large  displacements  and  changes  in  bonding  during  relaxation  leading  to 
a  substantial  lowering  of  the  surface  energy.  The  band  structure  is  also  obtained  for  the  lowest-energy  surface,  which  is  one  of  the  pos¬ 
sible  non-polar  terminations  of  the  (1  0  0).  The  results  provide  insight  into  the  growth  and  structure  of  f5-Ga203  nanoribbons. 
Published  by  Elsevier  B.V. 

Keywords:  Density  functional  theory;  Surface  structure;  Surface  energy;  Gallium  oxide;  Hartree-Fock;  Nanoribbons 


1.  Introduction 

Monoclinic  gallium  oxide  (P-Ga203)  has  a  number  of 
interesting  bulk  and  surface  properties.  When  doped 
(e.g.,  with  Sn)  it  can  be  used  to  form  an  ultraviolet- 
transmitting  (band  gap,  E^,  ^4.5^.9  eV)  conducting  layer 
[1].  Thin  films  of  polycrystalline  p-Ga203  containing  O 
vacancies  have  long  been  known  as  sensors  for  a  variety 
of  gases  including  H2  [2],  CH4  [3],  CO  [4]  and  O2  [5]  which 
change  the  electrical  conductivity  upon  adsorption.  A  fur¬ 
ther  source  of  interest  in  Ga203  surfaces  is  the  recent  work 
on  nanowires  and  nanoribbons  (or  nanobelts)  [6].  Nanorib¬ 
bons  grow  as  single  crystals,  typically  along  the  [0  1  0]  [7-9] 
or  [0  0  1]  [9-11]  directions;  although,  there  have  been 
reports  of  growth  along  the  [111]  [12],  [1  10]  [13]  or 
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[1  0  1]  [14]  directions.'  The  nanoribbons  have  well-defined 
surface  planes.  Growth  in  the  [0  0  1]  direction  gives  ribbons 
enclosed  by  ±(10  0)  and  ±(0  1  0)  planes,  while  ribbons 
grown  in  the  [0  1  0]  direction  have  sides  that  are  ±(  1  0  0) 
and  ±(10  1)  [6,9].  In  either  case,  the  wide  surfaces  of  the 
ribbon  are  ±(1  0  0).  Single  crystals  of  p-Ga203  cleave  read¬ 
ily  [15,16]  on  the  (1  0  0)  plane  which  suggests  that  this  is  the 
most  stable  surface. 

Furthermore,  powdered  Ga203  is  of  interest  in  catalysis, 
and  several  studies  have  used  the  infrared  absorption  [17- 
19]  or  electron  spin  resonance  [20]  spectra  of  adsorbates  on 
p-Ga203  powders,  activated  by  evacuation  at  elevated  tem¬ 
perature,  as  probes  of  the  surface  structure.  The  general 
conclusion  is  that  such  surfaces  exhibit  comparable  densi¬ 
ties  of  unsaturated  tetrahedral  and  octahedral  Ga  sites 
(see  below)  and  that  the  former  are  responsible  for  the 


*  Care  must  be  taken  in  reading  papers  referencing  the  |3-Ga203  crystal 
structure.  Some,  such  as  [14],  use  a  definition  of  the  unit-cell  axes  that 
differs  from  the  “standard”  convention  given  below. 
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Lewis-acid  character  of  the  activated  p-Ga203  powders.  In 
addition,  studies  of  the  thermal  oxidation,  in  dry  air  or  O2, 
of  hexagonal  GaN  epitaxial  films  [21,22]  powders  [22]  or 
nanowires  [23]  have  identified  P-Ga203  as  the  resulting 
surface  phase.  The  structure,  stability  and  chemical  proper¬ 
ties  of  the  oxide  are  significant  factors  in  the  subsequent 
processing  of  GaN  substrates  for  device  fabrication. 

In  view  of  this  rich  array  of  properties  and  applications 
it  appears  worthwhile  to  consider  in  detail  the  surface 
structure  of  P-Ga203.  Experimental  studies  of  the  growth 
and  ordering  of  well-characterized  P-Ga203  films  have 
been  reported  by  Franchy  and  co-workers  [24-30]  for  thin 
films  formed  by  oxidation  of  the  surfaces  of  the  crystalline 
bimetallic  compound  CoGa.  The  only  theoretical  studies  of 
P-Ga203  surfaces  of  which  we  are  aware  are  those  of  Kohl 
et  al.  [3 1  ]  and  Gonzalez  et  al.  [32]  (discussed  below)  which 
employed  semi-empirical  cluster  treatments  to  examine  the 
adsorption  of  CH4  and  H,  respectively,  on  (1  0  0)  surfaces. 
In  the  present  work,  ah  initio  calculations  will  be  done  to 
study  the  structure  and  charge  distribution  of  low-index 
surface  planes  of  P-Ga203  since,  as  is  already  known  [33], 
the  Ga203  surface  structure  can  have  a  pronounced  effect 
on  reactivity.  Attention  is  focused  on  the  (10  0),  (0  1  0), 
(0  0  1)  and  (1  0  I)  faces  since,  as  noted  above,  these  appear 
to  be  among  the  most  technologically  important.  There 
has,  to  our  knowledge,  been  no  previous  theoretical  or 
experimental  work  aimed  at  comparing  the  structure  and 
properties  of  these  low-index  P-Ga203  surface  planes. 

2.  Computational  details 

The  large  number  of  atoms  and  low  symmetry  of  the  unit 
cell  (see  below)  make  the  computation  demanding  in  terms 
of  machine  time.  Therefore,  geometry  optimizations  were 
performed  using  a  restricted  Hartree-Fock  (RHF)  ap¬ 
proach  with  Durand-Barthelat  large-core  pseudopotentials. 
The  electronic  properties  of  the  optimized  structures  were 
then  obtained  using  density  functional  theory  (DFT)  with 
all-electron  basis  sets.  Except  where  noted,  all  calculations 
were  done  using  the  CRYSTAL2003  suite  of  programs 
[34,35]  which  uses  Bloch  functions  constructed  from  local¬ 
ized  Gaussian  basis  sets.  Two-dimensionally-periodic  slabs, 
designed  to  maintain  the  symmetry  and  stoichiometry  of 
the  bulk  unit  cell,  were  used  to  model  the  various  surfaces. 
The  possible  role  of  adsorbates  and/or  defects  in  stabilizing 
polar  or  non-stoichiometric  surfaces  [36]  was  not  examined. 
Since  localized,  rather  than  plane-wave,  basis  sets  were  used 
there  was  no  need  to  impose  periodicity  along  the  surface 
normal. 

Optimization  of  the  basis  sets  and  the  unit-cell  lattice 
constants  was  done  using  the  LoptCG  script,^  which  com- 


^  The  LoptCG  script  was  obtained  at  http://www.crystal.unito.it/.  At 
present,  the  publicly-available  version  of  the  script  works  only  with 
CRYSTAL98.  The  latest  version,  CRYSTAL2003,  can  generate  gradients 
analytically  but  only  for  the  positions  of  individual  atoms,  not  for  basis-set 
parameters  or  lattice  constants. 


putes  gradients  numerically  using  CRYSTAL98,  followed 
by  the  application  of  a  conjugate-gradient  algorithm  to  lo¬ 
cate  the  total-energy  minimum.  Optimization  of  the  posi¬ 
tions  of  individual  atoms  was  done  using  CRYSTAL2003 
which  obtains  energy  gradients  analytically.  The  single¬ 
point  DFT  calculations  used  the  hybrid  B3LYP  functional 
which  has  been  shown  [37,38]  to  provide  more  accurate 
band  gaps  than  do  purely  ab  initio  functionals. 

The  pseudopotential  basis  sets  are  given  elsewhere.^  For 
the  Ga  21G*  set,  the  exponents  of  the  outer  sp  and  d  shells 
were  re-optimized  for  Ga203,  giving  0.2633  and 
0.2930  bohr“^  respectively.  For  the  O  41G  set,  re-optimiza- 
tion  of  the  outer  sp  exponent  gave  0.2178  bohr“^.  The  Ga 
all-electron  basis  set  was  the  86-41  lld41G  with  the  expo¬ 
nents  of  the  outer  three  sp  and  one  d  shells  re-optimized 
to  1.8605,  0.7105,  0.2775  and  0.6815  bohr“^,  respectively. 
The  O  all-electron  basis  set  was  the  8-41  IdlG  with  the 
exponents  of  the  outer  two  sp  and  one  d  shells  re-optimized 
to  0.4727,  0.1985  and  0.4497  bohr“^,  respectively. 

The  /f-point  sampling  used  an  8  x  8  x  6  grid  for  the  bulk 
and  8x6x  1,  6x2x  1,  8x4x  1  and  8x2x1  grids  for  the 
(1  0  0)-,  (0  1  0)-,  (0  0  1)-  and  (1  0  I)-(l  x  1)  slabs,  respec¬ 
tively.  In  treating  (2  x  2)  slab  supercells,  for  the  purpose 
of  including  longer-range  interactions,  the  k-point  grid 
was  reduced  by  a  factor  of  2  (e.g.,  a  4  x  3  x  1  grid  for  the 
(10  0)-(2x2)  slab).  The  interval  in  any  case  was  in  the 
range  of  0.10-0.14  bohr“'  per  point.  In  CRYSTAL,  trun¬ 
cation  of  the  sums  of  Coulomb  and  exchange  terms  in 
the  Fock  matrix  is  determined  by  five  overlap  criteria 
(T1-T5)  [34,35].  These  were  set  at  10“^  for  T1-T4  and 
10“'"^  for  T5.  For  integrals  involving  pseudopotentials  the 
truncation  criterion  was  typically  set  at  10“’.  In  the  multi¬ 
polar  expansion  zone  [34,35]  a  maximum  order  of  L  =  6 
was  used.  For  the  DFT  calculations  a  Fock  mixing  of 
50%,  a  Fermi-level  smearing  of  2  x  10“^  Hartrees  and  a  le¬ 
vel  shift  of  5  Hartrees  were  typically  used  in  order  to  aid 
self-consistent  field  (SCF)  convergence. 

3.  Results  and  discussion 

3.1.  Bulk  properties 

The  most  stable  polymorph  of  Ga203  is  the  monoclinic 
P  form  which  belongs  to  the  C2/m  (or  C2h^)  space  group 
(number  12)  with  lattice  parameters  and  atom  positions  gi¬ 
ven  in  Table  1  [39,40].  In  the  “standard”  setting,  which  is 
used  here,  the  h-axis  of  the  centrosymmetric  unit  cell  is 
the  twofold  rotational  axis,  and  the  mirror  plane  (i.e.,  the 
(0  1 0))  is  perpendicular  to  this  axis.  There  are  four 
Ga203  units  per  crystallographic  unit  cell  (which  comprises 
two  primitive  unit  cells)  with  two  inequivalent  Ga  sites  and 
three  inequivalent  O  sites.  Half  the  Gas  are  in  Ga(I)  sites 
which  form  slightly  distorted  tetrahedra  with  four  O  ions. 


^  The  basis  sets  (both  pseudopotential  and  all-electron)  with  full 
documentation  are  given  at  http://www.crystal.unito/it/Basis_Sets. 
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Table  1 

Observed  and  calculated  structural  parameters  for  P-Ga203 


Parameter 

Calculated^ 

Experimental 

a 

12.202  A 

12.23  ±0.02 

h 

3.035  A 

3.04  ±0.01 

c 

5.799  A 

5.80  ±0.01 

p 

103.70° 

103.7  ±0.3 

Ga(I)“ 

.v  =  0.08981;  z  =  -  0.20723 

.X  =  0.0904;  z  =  -  0.2052 

Ga(II) 

0.34195;  -0.31567 

0.3414;  -0.3143 

0(1) 

0.15754;  0.10867 

0.1674;  0.1011 

0(11) 

0.49553;  0.25810 

0.4957;  0.2553 

O(III) 

0.82544;  0.43290 

0.8279;  0.4365 

“  This  work,  obtained  in  an  RHF/pseudopotential  calculation  (see  text). 
The  crystallographic  unit-cell  volume  is  208.65  A^. 

Results  given  by  Geller  [39,  Table  II].  The  unit-cell  volume  is 
209.50  A^.  These  values  have  been  revised  slightly  by  Ahman  et  al.  [40]. 

“  The  atoms  are  all  at  the  4/  special  positions  of  the  C2/m  space  group, 
given  by  (0,0,0;  1/2, 1/2,0)  ±  (xOr).  For  each  inequivalent  site  the  values 
of  X  and  z  are  given. 

and  the  other  half  are  in  Ga(II)  sites  which  form  highly  dis¬ 
torted  octahedra  with  six  O  ions.  Each  0(1)  is  threefold 
coordinated  and  lies  at  the  intersection  of  two  octahedra 
and  one  tetrahedron.  Each  0(11)  is  also  threefold  coordi¬ 
nated  and  is  shared  between  one  octahedron  and  two  tetra- 
hedra.  Each  O(III)  is  fourfold  coordinated  and  lies  at  the 
corner  of  three  octahedra  and  one  tetrahedron.  Fig.  1 
shows  a  model  of  the  lattice,"*  which  illustrates  the  struc¬ 
tural  details  noted  above.  In  terms  of  the  crystallographic 
unit  cell  vectors  a,  b  and  c,  the  primitive  unit  cell  vectors 
are  given  by  a'  =  (a  —  b)/2,  b'  =  (a  +  b)/2  and  c'  =  c. 

Table  1  gives  the  lattice  constants  and  atom  positions 
obtained  in  the  present  work  using  the  RHF/pseudopoten¬ 
tial  structural  optimization  described  above.  The  optimiza¬ 
tion  was  done  iteratively  by  alternately  refining  the  lattice 
constants  and  the  atom  positions  in  separate  calculations 
until  convergence  was  achieved.  The  calculated  structure 
is  seen  to  be  in  good  agreement  with  experiment.  Table  2 
gives  the  various  Ga-O  nearest-neighbor  distances,  Mul- 
liken  overlap  populations  (OPs)  and  ionic  charges.  The 
Ga-O  distances  differ  from  experiment  [39]  by  <0.02  A. 
The  OPs  are  somewhat  smaller  than  those  found  in  Ex¬ 
tended  Hiickel  Theory  results  [41]  which  are  in  the  range 
of  0.26-0.29  for  Ga(I)-0  and  0.14-0.23  for  Ga(II)-0 
bonds.  However,  the  general  result  that  the  Ga(I)-0  bonds 
are  more  covalent  is  reproduced  in  the  present  calculations. 
The  ionic  charges  are  somewhat  larger  than  those  obtained 
[42]  in  tight-binding  calculations,  which  give  +1.80  and 
+1.64  for  Ga(I)  and  Ga(II),  respectively,  and  —1.09  to 
—  1.18  for  the  O  ions. 

The  bulk  band  structure  of  P-Ga203  has  been  studied 
theoretically  by  several  groups  [41^4].  The  bulk  valence 
band  density  of  states  was  computed,  using  DFT  with 
all-electron  basis  sets  as  described  above,  for  the  optimized 


Fig.  1  was  constructed  using  the  VENUS  set  of  programs  (R.A. 
Dilanian  and  F.  Izumi)  obtained  at  http://homepage.mac.com/fujioizumi/ 
vis  ualiza  tion/VENU  S  .html . 


Fig.  1.  The  [i-GazOs  lattice  viewed  along  the  (0  1  0)  direction.  The 
different  inequivalent  atom  types  are  labeled.  The  labels  “A”  and  “B” 
indicate  the  different  (10  0)  and  (0  0  1)  surface  terminations  (see  text),  and 
the  dashed  line  shows  the  (0  0  1)-B  plane.  The  box  shows  the  unit  cell,  and 
the  diagonal  plane  is  the  (1  0  I).  The  slightly  different  sizes  shown  for  the 
two  different  Gas  and  for  the  three  different  Os  are  not  significant. 

lattice  structure  given  in  Table  1.  The  results  (not  shown) 
are  in  good  agreement  with  previous  work^  [43]  which  used 
the  FLAPW  method  with  the  Perdew-Wang  exchange- 
correlation  functional.  The  band  structure  was  also  ob¬ 
tained  (not  shown)  and  found  to  be  in  good  agreement  with 
previous  work  [43].  In  particular  a  deep  conduction  band 
minimum  (CBM)  occurs  at  the  T-point,  and  the  valence 
band  maximum  (VBM)  is  nearly  flat.  The  band  gap  ob¬ 
tained  in  the  present  work  is  Eg  =  5.09  eV  vs.  experimental 
values  in  the  range  of  4. 5^.9  eV  [29,45,46].  A  calculated  Eg 
that  is  somewhat  larger  than  the  observed  value  is  typical 
[37]  for  the  B3LYP  functional. 

3.2.  The  (10  0)  surface 

Fig.  1  indicates  two  possible  (1  0  0)  surface  terminations, 
labeled  ”A”  and  “B”,  which  correspond  to  stoichiometric 
unit  cells  with  non-polar  surfaces.  These  are  described  as 
follows. 

A\  terminated  in  rows  of  0(II)’s  lying  along  the  [0  1  0] 
direction  with  each  0(11)  back-bonded  to  two  Ga(I)s. 
Between  each  0(11)  row  is  a  row  of  Ga(II)  sites.  The 
surface  0(11)  and  Ga(II)  atoms  are  each  singly- 
unsaturated,  i.e.,  missing  a  single  nearest  neighbor 
vs.  the  corresponding  bulk  sites.  Ga(I),  0(1)  and 
O(III)  atoms  at  the  surface  are  fully  coordinated,  as 
in  the  bulk.  This  surface  is  formed  by  cleaving  the 
Ga(II)-0(II)  bonds  shown  in  plane  A  in  Fig.  1,  while 
leaving  the  Ga(I)-0(II)  bonds  intact,  which  results  in 
two  identical  cleavage  surfaces. 


^  The  convention  used  in  [43]  for  labeling  the  O  ions  differs  from  that 
used  here  and  in  [39].  The  0(2)  [0(3)]  of  [43]  is  the  O(III)  [0(11)]  of  [39]. 
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Table  2 

Computed  bulk  nearest-neighbor  distances,  Mulliken  overlap  populations  and  ionic  charges^  for  p-Ga203 


0(1) 

0(11) 

O(III)’’ 

Charge 

Ga(I) 

1.822  (1)  [0.169] 

1.824  (2)  [0.162] 

1.850  (1)  [0.150] 

+2.080 

Ga(II) 

1.934  (2)  [0.118] 

1.932  (1)  [0.116] 

1.994  (1);  2.081  (2)  [0.094];  [0.080] 

+2.164 

Charge 

-1.412 

-1.383 

-1.449 

^  Entries  such  as  Ga(I)-0(I)  are  nearest-neighbor  distances  in  A.  Values  in  parentheses  are  the  numbers  of  nearest-neighbors  of  a  given  type,  and  charges 
are  in  units  of  \e\.  Numbers  in  brackets  are  Mulliken  overlap  populations  (OPs)  for  the  pair  of  nearest-neighbor  atoms.  For  example,  Ga(I)  has  two  0(11) 
nearest  neighbors  at  a  distance  of  1.824  A  with  a  bond  OP  of  0.162,  and  the  Ga(I)  ionic  charge  is  +2.080|e|. 

Each  Ga(II)  has  one  O(III)  at  1.994  A  and  two  more  at  2.081  A.  The  shorter  of  the  two  Ga(II)-0(III)  bond  lengths  corresponds  to  a  slightly  larger  OP. 


B:  terminated  in  nearest-neighbor  rows  of  Ga(II)  and 
O(III)  atoms,  each  singly-unsaturated.  Ga(I)  and 
0(1)  atoms  at  the  surface  are  fully  coordinated, 
and  there  are  no  0(11)  atoms  in  the  surface  plane. 

These  are  both  “Type  2”  structures  as  defined  by  Tasker 
[47,48]  in  that  each  individual  layer  is  charged,  but  the 
repeat  unit  (i.e.,  the  unit  cell)  has  no  static  dipole  moment. 
Hence,  the  surface  energy  (see  below)  is  expected  to  con¬ 
verge  to  a  finite  value  with  increasing  slab  thickness. 

The  (10  0)  surface  structures  used  in  previous  cluster 
calculations  differ  from  those  described  above.  Kohl  et  al. 
[31],  in  their  investigation  of  CH4  adsorption,  used  a  sur¬ 
face  like  A  but  with  a  complete  layer  of  0(11)  atoms  for 
which  all  Ga  sites  are  fully  coordinated.^  This  is  a  polar 
surface  since,  for  a  stoichiometric  lattice,  the  corresponding 
(Too)  surface  would  be  terminated  in  a  layer  composed  of 
rows  of  doubly-unsaturated  Ga(I)  and  singly-unsaturated 
Ga(II)  sites  (see  Fig.  1).  Gonzalez  et  al.  [32],  in  their  study 
of  H  adsorption,  in  fact  used  this  polar  Ga-terminated  sur¬ 
face  as  model  for  a  surface  with  O  vacancies.  In  the  absence 
of  charge-compensating  adsorbates  or  defects  [47,48]  this 
type  of  structure  would  not  be  expected  to  yield  a  finite  sur¬ 
face  energy. 

The  results  for  the  ( 1  0  0)-A  surface  will  be  discussed 
first.  The  calculations  were  done  using  slabs  of  A  =  20, 
40  and  60  layers,  where  A  =  20  corresponds  to  one  crys¬ 
tallographic  unit-cell  distance  in  the  [1  0  0]  direction  (cf. 
Fig.  1).  The  symmetry  and  stoichiometry  of  the  bulk  unit 
cell  are  thus  preserved,  and  there  is  only  one  atom  per 
layer  in  the  primitive  unit  cell.  In  view  of  the  results  of 
Gomes  et  al.  [49],  the  slabs  were  “cut”  from  the  opti¬ 
mized  bulk  lattice  described  above,  rather  than  from 
the  experimental  lattice.  Geometry  optimization  was 
done  globally,  i.e.,  by  allowing  complete  freedom  of  mo¬ 
tion  to  all  atoms  in  the  slab.  However,  the  symmetry  of 
the  unrelaxed  slab  was  maintained  so  that  displacements 
of  corresponding  atoms  in  layers  m  and  N  —  m+  \  were 
the  same.  Relaxation  of  the  whole  slab  (rather  than  of 
only  the  outermost  layer  or  two)  was  done  in  view  of 


^  The  notation  in  [31]  differs  from  that  used  here.  Surface  O  atoms  are 
labeled  “0(1)”  and  “0(11)”  according  to  the  type  of  Ga  site  to  which  they 
are  bonded.  However  (cf.  Fig.  1)  all  surface  O  atoms  in  this  model  are 
properly  termed  “0(11)”  in  the  present  usage. 


the  large  displacements,  extending  several  layers  into 
the  slab,  that  have  been  reported  for  Cr203  [50]  and 
for  various  forms  of  AI2O3  [51-54].  A  slab  is  considered 
sufficiently  thick  if  the  displacements  thus  generated  be¬ 
come  essentially  zero  in  all  directions  as  the  middle  of 
the  slab  is  approached. 

In  this  and  in  following  sections,  the  effects  of  relaxation 
will  be  described  in  three  complementary  ways.  First,  plots 
of  the  atomic  displacements  along  the  surface  normal  and 
in  the  surface  plane  will  be  given  vs.  depth  into  the  slab 
to  show  the  magnitude  of  the  effect  and  its  extent  into 
the  bulk.  Second,  the  displacements  will  be  shown  graphi¬ 
cally,  superimposed  on  the  ideally-terminated  unit  cell. 
However,  for  clarity,  only  the  outermost  layer  of  each 
inequivalent  atom  type  will  be  shown.  Third,  the  effects 
on  local  bonding  will  be  described  in  terms  of  ionic 
charges,  nearest-neighbor  distances  and  bond  OPs  relative 
to  the  corresponding  bulk  quantities. 

Fig.  2  shows  the  displacements  vs.  layer  number,  start¬ 
ing  from  the  surface  plane  (m=\)  and  progressing  into 
the  bulk,  for  the  optimized  A  =  60  slab.  The  displacements 
are  defined  with  respect  to  the  corresponding  positions  in 
the  bulk  lattice.  One  small  tick  mark  on  the  vertical  axis 
corresponds  to  0.02  A,  or  about  1%  of  a  typical  Ga-O 
nearest-neighbor  distance  (cf.  Table  2).  Thus,  the  effect  of 
relaxation  may  be  considered  to  be  small  for  m  ^  10 
(i.e.,  about  half-way  into  the  outermost  crystallographic 
unit  cell).  Fig.  3  shows  the  surface  relaxation  graphically. 
For  clarity,  the  displacement  vectors  have  been  magnified 
by  a  factor  of  4  relative  to  the  interatomic  distances. 

The  relaxation  process  is  seen  to  be  complex,  with  the 
outermost  Ga(I)  and  0(1)  atoms  moving  up  toward  the 
surface  and  the  0(11)  surface  layer  displacing  laterally. 
However,  there  is  little  difference  between  the  bulk  and 
surface  ionic  charges  which,  for  corresponding  atom 
types,  differ  by  ^0.06|e|  (cf.  Table  2).  As  expected,  the 
greatest  difference  between  bulk  and  surface  ion  charges 
occurs  for  the  unsaturated  Ga(II)  and  0(11)  sites  which 
are  slightly  less  ionic  than  in  the  bulk.  Relative  to  the 
bulk  (Table  2),  the  Ga(I)-0(I)  and  Ga(I)-0(III)  distances 
are  longer  (1.911  and  2.027  A,  respectively)  and  the 
Ga(I)-0(II)  distance  shorter  (1.720  A).  The  Ga(II)-0(I) 
distance  is  virtually  unchanged  (1.932  A)  while  the 
Ga(ir)-0(III)  distances  are  slightly  shorter  (1.946  and 
2.024  A).  Thus,  the  Ga(I)  tetrahedron  appears  to  be  more 
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Fig.  2.  Displacements  (in  A),  vs.  layer  number,  along  the  surface  normal 
[(5(z)]  and  in  the  surface  plane  [6(x,y)'\  relative  to  the  ideally-terminated 
positions  for  the  { 1  0  0)-A  surface.  The  results  are  those  computed  for  a 
fully-relaxed  N  =60  slab.  Layer  number  1  is  the  outermost  surface  layer, 
and  S(z)  >  0  corresponds  to  an  outward  displacement. 


[001] 


a  (100)-4  b 

Layer  # 


Fig.  3.  Schematic  model  of  the  ideally-terminated  (1  0  0)-A  surface,  (a) 
The  view  along  the  surface  normal  with  the  dashed  lines  indicating  the 
surface  unit  cell.  For  simplicity,  only  the  outermost  layer  of  each  type 
(0(1),  Ga(I),  etc.)  is  shown,  (b)  The  outermost  surface  layers,  viewed  along 
the  [0  1  0]  direction,  with  vectors  indicating  the  displacements  occurring 
during  relaxation.  For  clarity,  the  displacement  vectors  have  been 
magnified  by  a  factor  of  4  relative  to  the  interatomic  distances.  All 
displacements  are  in  the  plane  of  the  page,  i.e.,  the  (0  1  0)  plane.  For 
simplicity,  displacements  of  translationally-equivalent  atoms,  which  are 
identical  to  those  shown,  have  been  omitted. 


distorted  on  the  relaxed  surface  than  in  the  bulk.  Qualita¬ 
tively,  the  OPs  increase,  vs.  the  bulk  values,  as  the  bond 
distances  decrease. 

The  surface  energy  of  a  two-dimensionally-periodic  slab 
with  a  thickness  of  K  primitive  unit  cells,  Gk,  is  defined  as 
[51]  E(K) —  2A(Tic+ KE^,  where  E(K)  is  the  total  energy 
per  slab  unit  cell,  E^  is  the  energy  per  unit  cell  of  the  bulk 
and  A  is  the  area  of  the  slab  unit  cell.  Both  faces  of  the  slab 
are  taken  to  be  identical.  The  practical  difficulties  in 
obtaining  a  converged  value,  o  —  oo),  and  in 

choosing  a  self-consistent  E^  are  discussed  elsewhere  [55]. 
In  the  present  work,  E{K)  was  computed  in  single-point 
DFT  calculations  for  both  ideal  and  globally-relaxed  slabs 
(see  above)  with  N  =  20,  40  and  60  (i.e.,  Ai  =  2,  4  or  6  prim¬ 
itive  unit  cells)  and  found  to  be  linear  vs.  K  over  this  range, 
which  implies  that  convergence  is  achieved  for  N  ^  20.  A 
value  for  a  (either  ideally-terminated  or  relaxed)  was  then 
obtained  from  the  zero-intercept  of  the  linear  least-squares 
fit,  and  the  results  are  given  in  Table  3. 

The  calculations  described  thus  far  are  based  on  a  ( 1  x  1) 
surface  unit  cell.  To  determine  whether  the  inclusion  of 
longer-range  interactions  would  yield  a  further  lowering 
of  energy,  a  calculation  was  done  for  an  iV  =  20  slab  with 
a  (2  X  2)  unit  cell.  This  allows  for  the  possibility  of  recon¬ 
struction  as  well  as  relaxation.  However,  the  optimized 
RHF/pseudopotential  total  energy  of  the  fully-relaxed 
(2x2)  structure  differed  from  that  of  the  (1x1)  by  less 
than  1  meV  per  (1x1)  slab  unit  cell.  Hence,  the  surface 
structure  is  adequately  modeled  with  the  (1x1)  unit  cell. 

Similar  calculations  were  done  for  a  ( 1  x  1 )  unit  cell  with 
the  (1  0  0)-B  termination,  and  the  results  are  shown  in  Figs. 
4  and  5.  The  surface  energy  is  even  lower  for  this  structure 
than  for  the  ( 1  0  0)-A.  The  displacements  of  atoms  for  the 
(1  0  0)-B  are  generally  smaller  in  magnitude  than  are  those 
of  corresponding  atoms  for  the  ( 1  0  0)-A,  which  is  consis¬ 
tent  with  (1  0  0)-B  being  the  more  stable  of  the  two  before 
relaxation.  Also,  the  3{z)  values  indicate  a  predominantly 
inward  displacement  of  atoms  near  the  (1  0  0)-B  surface. 
The  fact  that  the  relaxed  ( 1  0  0)-B  surface  has  the  lowest 
energy  of  those  considered  here  (cf.  Table  3)  is  consistent 
with  the  observation  [15,16]  that  P-Ga203  cleaves  easily 


Table  3 

Surface  energies  {a,  J  m“^)  for  the  ideally-terminated  and  the  relaxed  or 
reconstructed  surfaces 


Area 

Ideal  a 

Relaxed  a 

(1  0  0)-A 

0.1760 

1.68 

1.13^* 

(1  0  0)-B 

0.96 

0.68‘‘ 

(0  10) 

0.6874 

2.78 

2.03^ 

(0  0  1)-A 

0.1852 

3.35 

(0  0  1)-B 

2.65 

1.40’’ 

(10  1) 

0.3705 

3.99 

1.57'" 

The  area  (nm^)  is  that  of  the  primitive  (1x1)  surface  unit  cell.  The  labels 
“A”  and  “B”  refer  to  different  lattice  terminations. 

“  Relaxed  (1x1)  unit  cell.  1  eV  nm“^  =  0.160219  J  m“^. 

’’  Reconstructed  (2  x  2)  unit  cell.  A  minimum-energy  configuration  could 
not  be  located  by  relaxing  a  (1  x  1)  unit  cell  (see  text). 
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on  the  (10  0)  plane.  The  (1  0  0)-B  surface^  is  also  the  one 
proposed  [24,28,29]  for  crystalline  P-Ga203  films  formed 
by  oxidation  of  CoGa  (10  0)  and  (0  0  1)  surfaces.  It  has 
also  been  noted  [9]  that  P-Ga203  nanoribbons  grow  with 
( 1  0  0)  as  the  broad  surface  which  is  consistent  with  the  rel¬ 
atively  low  a  for  the  (1  0  0)-B.  The  greater  stability  of  the 
( 1  0  0)-B  vs.  -A  surface  is  probably  due  to  the  lesser  degree 
of  unsaturation  for  the  former.  On  ( 1  0  0)-B  the  surface 
O(III)  atoms  are  missing  one  out  of  four  Ga  nearest-neigh- 
bors;  whereas,  on  ( 1  0  0)-A  the  surface  0(11)  atoms  are 
missing  one  out  of  three  bulk  nearest-neighbor  Ga  atoms. 

In  terms  of  local  bonding,  the  relaxed  (1  0  0)-B  surface 
differs  only  slightly  from  the  bulk  which  is  consistent  with 
its  relatively  low  a.  The  surface  ionic  charges  differ  from 
those  in  the  bulk  by  <0.05|e|.  The  Ga(I)-0(I)  and  Ga(I)- 
0(11)  distances  (1.843  and  1.839  A,  respectively)  are  slightly 
larger  than  in  the  bulk,  and  the  Ga(I)-0(III)  distance  is 


’  The  definition  of  the  crystal  axes  used  in  [24,28,29]  differs  from  that  in 
[39,40]  and  in  the  present  work.  The  axes  a,  h,  c  in  the  former  correspond 
to  the  axes  b,  c,  a  in  the  latter.  Hence,  the  (0  0  1)  plane  described  in 
[24,28,29]  is  the  (10  0)  plane  of  the  present  work. 


Fig.  5.  Similar  to  Fig.  3  but  showing  results  for  the  (1  0  0)-B  surface.  In 
(a),  only  Ga(II)  and  O(III)  exhibit  a  significantly  large  relaxation 
displacement.  Other  displacements  are  too  small  to  be  seen  clearly,  even 
when  magnified,  and  have  been  omitted  from  (b). 

slightly  smaller  (1.844  A),  indicating  a  more  regular  tetra¬ 
hedron  at  the  surface.  The  Ga(II)-0(I)  distance  (1.978  A) 
is  slightly  greater  than  in  the  bulk,  and  the  Ga(II)-0(lI) 
distance  is  slightly  smaller  (1.887  A).  The  distance  from 
Ga(II)  to  the  two  equidistant  O(III)  sites  is  1.974  A,  slightly 
less  than  the  bulk  value  of  2.081  A. 

3.3.  The  (010)  surface 

For  this  surface  there  is  only  one  termination,  which 
contains  all  five  atom  types  with  two  Ga203  units  per  sur¬ 
face  unit  cell  per  slab  layer.  Thus,  the  (0  1  0)  orientation  in¬ 
volves  10  atoms  per  layer  unlike  the  (10  0),  discussed 
above,  with  only  one  atom  per  layer.  Figs.  6  and  7  give 
the  structural  results  for  this  surface,  which  can  be  de¬ 
scribed  as  follows: 

[Ga( I)-.]  bonded  to  one  0(1)  and  one  O(III)  in  plane;  one 

0(11)  below. 

[Ga)!!)-.]  bonded  to  one  0(11)  and  one  O(III)  in  plane; 

one  0(1)  and  one  O(III)  below. 

[0(!y.]  bonded  to  one  Ga(I)  in  plane;  one  Ga(II)  below. 

[0(!!y.]  bonded  to  one  Ga(Il)  in  plane;  one  Ga(I)  below. 

[0(!!!y.]  bonded  to  one  Ga(I)  and  one  Ga(II)  in  plane; 

one  Ga(II)  below. 

All  surface  atoms  are  singly-unsaturated  except  for  the 
Ga(II)  which  is  doubly-unsaturated.  The  surface  layer 
may  be  viewed  (Fig.  7(a))  as  consisting  of  individual 
Ga203  units  back-bonded  to  underlayer  atoms,  with  no 
direct  interaction  between  surface  Ga203  units. 

Global  geometry  optimizations  were  performed  for 
A  =  5,  7,  1 1  and  15  with  convergence  being  achieved,  as  de¬ 
fined  above,  for  N  ^  1.  The  displacements  obtained  for 
A  =  15  are  given  in  Figs.  6  and  7(b),  and  the  surface  ener¬ 
gies  for  the  ideal  and  the  relaxed  slabs  are  listed  in  Table  3. 
The  effect  of  relaxation  is  seen  to  be  small,  as  defined 
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Fig.  6.  Similar  to  Fig.  2  but  showing  results  for  the  (0  1  0)  surface  for  an 
=  15  slab.  The  “(5{x,  j')”  quantity  plotted  in  this  case  is  ((5{x)^  + 

Lines  connecting  the  O  points  are  visual  aids. 


(010) 


Fig.  7.  Similar  to  Fig.  3  but  showing  the  (0  1  0)  surface.  In  (a),  heavy  lines 
connecting  a  few  of  the  atoms  illustrate  the  surface  bonding  described  in 
the  text.  In  (b),  all  atoms  shown  are  in  the  same  (outermost)  plane  for  the 
ideally-terminated  surface.  The  0(1)  displacements  are  too  small  to  be  seen 
clearly. 


above,  beyond  the  third  or  fourth  layer  into  the  bulk  of  the 
slab,  where  three  layers  correspond  to  one  crystallographic 
unit  cell  distance  in  the  [0  1  0]  direction.  The  surface  ionic 


charges  were  found  to  be  somewhat  smaller  than  in  the 
bulk,  indicating  a  higher  covalency  at  the  surface.  The  re¬ 
laxed  N  =  \5  (0  10)  surface  charges  were  (cf.  Table  2) 
+2.017|e|  and  +2.09 l|e|  for  Ga(I)  and  Ga(II),  respectively, 
and  -1.358|e|,  -1.336H  and  -1.377|e|  for  0(1),  0(11)  and 
O(III),  respectively.  At  the  surface,  the  distances  from 
Ga(I)  to  0(1),  0(11)  and  O(III)  (1.707,  1.803  and  1.762  A, 
respectively)  are  all  shorter  than  in  the  bulk,  with  corre¬ 
spondingly  larger  OPs.  Likewise,  the  distances  from  Ga(II) 
to  0(1)  and  0(11)  (1.879  and  1.754  A)  are  shorter,  with  lar¬ 
ger  OPs,  than  in  the  bulk.  At  the  surface,  there  are  two  un¬ 
equal  Ga(II)-0(III)  bonds  which  are  again  shorter  (1.847 
and  1.910  A)  than  in  the  bulk.  The  results  indicate  an  in¬ 
creased  covalent  interaction  between  pairs  of  surface  atoms 
so  that  the  (0  1  0)  surface  resembles  to  some  extent  a  layer 
of  Ga203  “molecules”  chemisorbed  on  the  surface  of  the 
bulk  lattice. 

As  in  the  case  of  the  ( 1  0  0)-A  surface  discussed  above,  a 
geometry  optimization  was  also  done  for  an  A  =  5  slab 
with  a  (2x2)  supercell  (containing  a  total  of  200  atoms 
per  slab  unit  cell)  in  order  to  determine  if  reconstruction 
would  further  lower  the  total  energy  of  the  relaxed  struc¬ 
ture.  However,  the  RHF/pseudopotential  energy  of  the 
fully-relaxed  (2  x  2)  structure  differed  from  that  of  the 
(1  X  1)  by  less  than  1  meV  per  (1x1)  slab  unit  cell. 

3.4.  The  (001)  surface 

As  in  the  case  of  the  (1  0  0),  two  non-polar  (0  0  1)  termi¬ 
nations  of  the  stoichiometric  unit  cell  are  possible. 

A',  terminated  in  rows  of  singly-unsaturated  0(1)  atoms 
lying  along  the  [0  1  0]  direction  and  back-bonded  to 
two  Ga(II)s.  The  0(1)  rows  are  separated  by  rows 
of  singly-unsaturated  Ga(I)s.  Ga(II)  and  0(11)  atoms 
in  the  surface  plane  are  fully  coordinated,  and  there 
are  no  0(III)s  in  the  plane. 

B\  terminated  in  rows  of  doubly-unsaturated  O(III) 
atoms  back-bonded  to  one  Ga(I)  and  one  Ga(n) 
and  separated  by  rows  of  doubly-unsaturated  Ga(II) 
sites.  Ga(I)  and  0(11)  atoms  in  the  surface  plane  are 
fully  coordinated,  and  there  are  no  0(I)s  in  the  sur¬ 
face  plane. 

When  relaxing  (1x1)  unit  cells  for  either  termination  a 
minimum-energy  configuration,  in  which  the  forces  on  all 
atoms  are  effectively  zero,  could  not  be  located.  The  behav¬ 
ior  during  optimization  suggests  that  the  potential  energy 
surface  is  relatively  flat  near  the  energy  minimum,  and 
increasing  the  precision  limits  for  SCF  convergence  and 
for  the  truncation  of  the  Coulomb  and  exchange  terms 
(see  above)  did  not  rectify  the  situation.  However,  conver¬ 
gence  was  achieved  by  increasing  the  slab  unit  cell  from 
(1x1)  to  (2x2).  In  view  of  the  large  number  of  atoms 
involved  (e.g.,  160  atoms  for  a  (2x2)  A=  40  slab),  only 
the  B  termination  was  considered  since  a  for  the  ideally- 
terminated  (0  0  1)-B  is  significantly  lower  than  that  of  the 
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(0  0  1)-A  (Table  3).  Thus  it  is  assumed  that  the  (0  0  1)-B 
will  remain  the  lower  of  the  two  in  energy  even  for  a 
(2x2)  surface  unit  cell.  The  results  for  the  (OOl)-B 
(2  X  2)  surface  are  shown  in  Figs.  8  and  9.  The  distinction 
between  the  (2x2)  and  (1x1)  surface  unit  cells  is  very 
small,  in  that  there  are  only  slight  differences  (<0.0  0  1  A) 
in  the  displacements  of  atoms  in  the  (2  x  2)  cell  that  are 
translationally  equivalent  in  the  (1x1).  Hence,  there  is  lit¬ 
tle  or  no  indication  of  an  actual  reconstruction,  and  Fig.  9 
shows  only  the  (1x1)  cell. 

The  charges  on  the  doubly-unsaturated  surface  Ga(II) 
and  O(III)  ions  (+2.063|e|  and  — 1.317|e|,  respectively)  are 
smaller  in  magnitude  by  ^0.1  |e|  than  on  the  corresponding 
bulk  ions  (cf.  Table  2),  indicating  a  higher  degree  of  cova¬ 
lency.  However,  the  charges  on  the  other  surface  ions  differ 
by  at  most  0.016|e|  from  the  bulk  counterparts.  The 
Ga(I)-0(I)  and  Ga(I)-0(III)  distances  (1.799  and 
1.765  A,  respectively)  are  shorter  than  in  the  bulk;  whereas, 
the  Ga(I)-0(II)  distance  is  longer  (1.885  A).  Thus,  the 
Ga(I)  tetrahedron  is  more  distorted  than  in  the  bulk.  The 
Ga(II)-0(I)  and  Ga(II)-0(II)  distances  (1.837  and 
1.900  A,  respectively)  are  both  shorter  than  in  the  bulk. 


5  to  15  20 


Fig.  8.  Similar  to  Fig.  2  but  showing  results  for  the  (OOl)-B  (2x2) 
surface  for  an  =  40  slab.  As  noted  in  the  text,  the  (2  x  2)  surface  unit  cell 
was  used  as  a  means  of  promoting  convergence  of  the  geometry 
optimization,  but  there  is  little  or  no  indication  of  an  actual 
reconstruction. 


(OOl)-fl 

a  b 


Fig.  9.  Similar  to  Fig.  3  but  showing  results  for  the  (0  0  1)-B  surface.  Only 
the  (1x1)  primitive  and  crystallographic  surface  unit  cells  are  shown  (see 
text).  The  Ga(II)  and  O(III)  surface  sites  are  each  doubly-unsaturated  (i.e., 
missing  two  of  the  bulk  nearest-neighbors).  All  other  species  are  fully 
coordinated. 

The  one  Ga(II)-0(III)  bond  is  much  shorter  than  in  the 
bulk  (1.748  vs.  1.994  A)  with  a  much  higher  OP  (0.205 
vs.  0.094,  cf.  Table  1).  Thus,  the  Ga(II)  is  more  covalently 
bonded  at  the  surface  than  in  the  bulk. 

3.5.  The  {101)  surface 

The  (10  1)  surface  was  also  analyzed  in  view  of  its 
importance,  noted  above,  in  the  growth  of  P-Ga203  nano¬ 
ribbons.  This  surface  consists  of  (cf.  Fig.  1)  rows  of  doubly- 
unsaturated  0(1)  and  singly-unsaturated  O(III)  atoms  lying 
along  the  [0  1  0]  direction.  A  row  of  singly-unsaturated 
Ga(I)  sites  lies  adjacent  to  the  0(1)  row,  and  a  row  of  dou¬ 
bly-unsaturated  Ga(II)  sites  lies  adjacent  to  the  O(III)  row. 
The  surface  0(11)  atoms  are  fully  coordinated.  In  forming 
the  surface,  Ga(I)  loses  its  O(III)  bulk  nearest-neighbor, 
and  Ga(II)  loses  its  two  0(1)  bulk  nearest-neighbors  (cf. 
Fig.  1). 

Global  geometry  optimizations  were  done  for  N  =  20, 
40,  60  and  80,  with  convergence  being  achieved  for 
N  ^  40.  The  displacements  obtained  for  A  =  80  are  shown 
in  Figs.  10  and  11.  In  this  case,  relaxation  results  in  larger 
displacements  near  the  surface  than  were  found  for  other 
planes;  however,  these  become  negligibly  small  beyond 
m  «  20.  The  surface  energy  of  the  ideally- terminated  sur¬ 
face  (Table  3)  is  relatively  high,  as  expected  from  the  high 
degree  of  unsaturation  noted  above.  Relaxation  yields  a 
substantial  reduction  in  <7,  by  a  factor  of  ^2.5,  which  is  a 
larger  effect  than  for  any  other  surface  studied  here.  The 
relaxation  is  also  accompanied  by  a  larger  redistribution 
of  ionic  charges  than  was  noted  for  other  surfaces.  Table 
4  shows  the  nearest-neighbor  distances,  OPs  and  ionic 
charges  for  the  (1  0  I)  surface  in  the  same  format  as  those 
for  the  bulk  given  in  Table  2.  As  was  done  for  other  sur¬ 
faces,  a  calculation  was  performed  for  a  (2  x  2)  surface  unit 
cell  (with  N  =  40),  and  a  negligible  effect  on  the  total  en¬ 
ergy  was  again  found. 

The  origin  of  the  remarkable  stabilization  of  the  (1  0  I) 
surface  can  be  seen  by  comparing  Tables  2  and  4  and  by 
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Fig.  10.  Similar  to  Fig.  2  but  showing  results  for  the  {1  0  I)  surface  for  an 
=  80  slab.  Note  the  much  larger  magnitudes  for  displacements  near  the 
surface  in  comparison  to  Figs.  2,  4,  6  and  8. 


examining  Fig.  1 1 .  Fig.  1 1  is  only  indicative,  however,  since 
underlayer  atoms,  which  also  undergo  significant  displace¬ 
ments  (cf.  Fig.  10),  have  been  omitted  for  clarity.  As  noted 
above,  Ga(I)  and  0(1)  are  singly-  and  doubly-unsaturated, 
respectively,  on  the  ideally-terminated  surface.  However, 
relaxation  results  in  displacement  of  0(1)  and  Ga(I)  to¬ 
wards  each  other  to  form  a  “bulk-like”  bond 
(OP  =  0.101),  leaving  the  0(1)  fully  coordinated  and  the 
Ga(I)  fivefold  coordinated  (vs.  fourfold  in  the  bulk).  This 
in  turn  leads  to  a  weakening  of  the  Ga(r)-0(II)  bond 
(i.e.,  a  decrease  in  OP  from  0.162  in  the  bulk  to  0.125  at 
the  surface)  and  a  strengthening  of  the  Ga(II)-0(II)  inter¬ 
action  (i.e.,  an  increase  in  OP  from  0.116  to  0.154).  The 
Ga(II)-0(III)  interaction  at  the  surface  is  also  noticeably 
stronger  than  in  the  bulk. 
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©-©  o 
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Fig.  11.  Similar  to  Fig.  3  but  showing  results  for  the  (1  0  I)  surface.  The 
displacement  vectors  have  been  magnified  by  a  factor  of  2  relative  to  the 
interatomic  distances.  Note  that  in  corresponding  figures  for  other 
surfaces  the  magnification  is  by  a  factor  of  4.  The  O(III)  displacement  is 
too  small  to  be  seen  clearly. 


3.6.  The  (10  0)  surface  electronic  structure 

As  noted  previously,  bulk  single  crystals  of  p-Ga203 
cleave  easily  on  the  (10  0)  plane  [15,16].  Hence,  it  should 
be  possible  to  observe  experimentally  the  electronic  struc¬ 
ture  of  this  surface  by  cleaving  a  crystal  in  ultra-high  vac¬ 
uum  (UHV)  and  obtaining  ultraviolet  photoemission 
spectroscopy  (UPS)  data.  Charging  of  the  insulating  sam¬ 
ple  during  UPS  could  be  eliminated  by  first  annealing  in 
UHV  to  create  oxygen  vacancies,  thus  inducing  n-type 
semiconductor  behavior  [56].  A  UPS  experiment  of  this 
nature  has  not,  to  our  knowledge,  been  reported  yet.  It  is 
also  noted  that  the  ( 1  0  0)-B  structure  has  been  proposed 
for  crystalline  P-Ga203  films  grown  by  oxidation  of 
CoGa(l  0  0)  and  (0  0  1)  surfaces  [24,29];  although,  again, 
no  UPS  data  have  been  reported. 

In  order  to  examine  the  possibility  of  intrinsic  surface 
states,  the  surface  band  structure  was  computed  for  the 
fully-relaxed  (1x1)  N  —6D  slab  with  the  (1  0  0)-B  struc¬ 
ture,  and  the  results  (near  the  band  edges)  are  shown  in 
Fig.  12.  The  surface  Brillouin  zone,  which  is  simple  for 
the  (10  0)  plane,  was  constructed  using  the  monoclinie  re- 
eiprocal  lattiee  vectors  given  by  Hunderi  [57].  The  results 
show  basically  the  same  nearly-fiat  VBM  and  deep  CBM 
at  r  as  is  seen  for  the  bulk  [43],  and  there  is  no  indication 
of  surface-related  features  in  the  band  gap.  The  absence  of 
intrinsic  surface  states  near  the  VBM  should  assist  in  the 


Table  4 

Similar  to  Table  2  but  showing  results  for  the  relaxed  {1  0  1)  surface 


0(1) 

0(11) 

O(III) 

Charge 

Ga(I) 

Ga(II) 

Charge 

1.852  (1);  2.067  (2)  [0.189];  [0.101] 

-1.387 

1.989  (2)  [0.125] 

1.938  (1)  [0.154] 
-1.391 

1.848  (1);  2.072  (2)  [0.178];  [0.110] 
-1.346 

-f  2.1 18 

-f  2.056 

See  text  for  details. 
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Fig.  12.  Surface  band  structure  diagram,  computed  using  DFT  with  all¬ 
electron  basis  sets,  for  the  fully-relaxed  (1  x  1)  A'^  =  60  slab  with  the  (10  0)- 
B  termination.  The  inset  shows  the  surface  Brillouin  zone.  The  region 
shown  comprises  the  upper  (lower)  part  of  the  valence  (conduction)  band, 
and  the  zero  of  energy  is  at  the  valence  band  maximum.  The  y-axis  lies 
parallel  to  the  [0  1  0]  direction. 


detection  and  identification  of  extrinsic  surface  states  aris¬ 
ing  from  defects.  In  their  electron  energy  loss  data  for  P- 
Ga203  (1  0  0)  films  on  CoGa,  Schmitz  et  al.  [29]  reported 
a  peak  of  uncertain  origin  at  3.3  eV  in  addition  to  the  inter¬ 
band  transition  at  about  4.5  eV.  The  present  results  indi¬ 
cate  that  the  3.3  eV  feature  is  not  associated  with 
intrinsic  surface  states. 

4.  Summary 

In  conclusion,  the  physical  and  electronic  structure  of 
the  (1  0  0),  (0  1  0),  (0  0  1)  and  (1  0  I)  faces  of  P-Ga203  have 
been  addressed  using  ab  initio  theory.  For  the  bulk,  good 
agreement  is  obtained  with  available  experimental  data 
(i.e.,  crystal  structure  and  band  gap)  and  with  previous 
ab  initio  studies  (i.e.,  band  structure).  The  results  provide 
a  physical  picture  of  the  structure  of  each  surface  and  of 
how  that  structure  is  derived  by  relaxing  the  corresponding 
ideally-terminated  surface.  In  all  cases,  especially  for  the 
(1  0  I)  plane,  large  displacements  from  the  ideal  positions 
are  found  for  atoms  near  the  surface;  however,  such  effects 
decay  quickly  with  distance  into  the  bulk. 

The  stability  of  various  surface  planes  and  terminations 
has  been  assessed  via  calculation  of  the  respective  surface 
energies.  The  relatively  low  c  of  the  ( 1  0  0)-B  surface  is  con¬ 
sistent  with  the  observed  tendency  of  p-Ga203  single  crys¬ 
tals  to  cleave  on  the  (10  0)  plane  and  with  the  tendency  of 
P-Ga203  nanoribbons  to  grow  with  the  (10  0)  plane  as  the 
wide  surface.  Growth  with  this  structure  allows  nanorib¬ 


bon  formation  to  proceed  with  a  minimal  increase  in  sur¬ 
face  energy.  The  stability  of  the  (lOO)-B  surface  is 
ascribed  to  a  local  bonding  configuration  that  is  not  very 
different  from  that  of  the  bulk. 

For  growth  in  the  [0  0  1]  direction,  the  narrow  face  of 
the  nanoribbon  is  the  more  energetic  (0  1  0)  plane.  The  lo¬ 
cal  bonding  on  the  (0  1  0)  surface  is  found  to  resemble 
somewhat  that  of  a  layer  of  chemisorbed  Ga203  “mole¬ 
cules”.  However,  growth  in  the  [0  1  0]  direction  occurs  with 
(1  0  1),  and  not  (0  0  1),  as  the  narrow  face.  This  is  at  first 
surprising,  in  view  of  the  high  cr  of  the  ideally-terminated 
(1  0  I)  surface.  However,  Table  3  shows  that  relaxation 
leads  to  a  very  substantial  reduction  in  <7  for  the  (1  0  1), 
to  a  value  only  ^12%  larger  than  that  of  the  (0  0  1)-B. 
The  stabilization  of  the  (1  0  I)  is  analyzed  in  terms  of  a  sig¬ 
nificant  reconfiguration  of  the  local  bonding  at  the  relaxed 
surface.  It  is  possible  that  some  lesser  effect  not  considered 
here,  such  as  defect  formation,  further  stabilizes  the  (1  0  I) 
to  a  point  where  it  is  energetically  favored  over  the  (0  0  1) 
as  the  narrow  ribbon  face. 

Of  the  surfaces  considered,  the  most  stable  structures 
(i.e.,  (1  0  0)-B,  (1  0  0)-A  and  (0  0  1)-B,  in  order  of  increas¬ 
ing  a)  all  have  unsaturated  Ga(II)  sites  but  no  unsaturated 
Ga(I)  sites.  This  is  in  apparent  conflict  with  experimental 
results  [17]  for  Ga203  powders  which  indicate  essentially 
equivalent  bulk  and  surface  compositions.  For  p-Ga203 
this  means  that  the  powder  surface  has,  on  average,  nearly 
equal  concentrations  of  unsaturated  Ga(I)  and  Ga(ir)  sites. 
One  factor  could  be  the  tendency  to  form  O  vacancies  dur¬ 
ing  heating  in  vacuo  [56],  which  is  done  to  prepare  clean 
Ga203  powder  surfaces  [17].  However,  an  impossibly  high 
concentration  of  such  vacancies  would  be  needed  to  ac¬ 
count,  by  this  means  alone,  for  the  experimental  observa¬ 
tion.  Another  possibility  is  that  the  (1  0  I)  surface,  which 
has  equal  concentrations  of  unsaturated  Ga(I)  and  Ga(II), 
might  be  favored  for  p-Ga203  in  the  form  of  small  particles 
as  it  is  for  nanoribbons  grown  in  the  [0  1  0]  direction.  A 
further  possibility  is  that,  for  example,  the  polar  (10  0)  sur¬ 
face  termination  consisting  of  a  layer  of  unsaturated  Ga(I) 
and  Ga(II)  sites  [32]  might  be  stabilized  on  powder  surfaces 
by  the  formation  of  charge-compensating  defects  [48]. 
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